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ABSTRACT 

We report on the results of a simultaneous monitoring campaign employing eight Chandra X-ray 
(0.5-10 keV) and six VLA/EVLA (8.4 GHz) radio observations of NGC 4051 over seven months. 
Evidence for compact jets is observed in the 8.4 GHz radio band; This builds on mounting evidence 
that jet production may be prevalent even in radio-quiet Seyferts. Assuming comparatively negligible 
local diffuse emission in the nucleus, the results also demonstrate an inverse correlation of -Lradio c< 
-^x-ray^'' ^ Configuration is excluded in the case that diffuse emission plays a significant 

role, the relation is still £radio « ix-ray^^"^- Current research linking the mass of supermassive 
black holes and stellar-mass black holes in the "low/hard" state to X-ray luminosities and radio 
luminosities suggest a "fundamental plane of accretion onto black holes" that has a positive correlation 
of iradio cx ^x-^ay ^ur simultancous results differ from this relation by more than II a {6 a 
excluding the A configuration), indicating that a separate mode of accretion and ejection may operate 
in this system. A review of the literature shows that the inverse correlation seen in NGC 4051 is seen 
in three other black hole systems, all of which accrete at near 10% of their Eddington luminosity, 
perhaps suggesting a distinct mode of disk-jet coupling at high Eddington fractions. We discuss our 
results in the context of disk and jets in black holes and accretion across the black hole mass scale. 



1. INTRODUCTION 

With evidence of supermassive black holes 
(SM BH) lurking at the center of nearly all galax- 
ies (Richs tone et al.l |1998[ ). it is pertinent to examine 
and understand their properties as well as their impacts. 
Studies have shown a critical relation between SMBHs 
and their host galaxies in the form of the M-L and 
M-a relations. The M-L relation between the mass of 
the SMBH and the luminosity of the bulge suggests an 
intrinsic link between the SMBH and the amount of mass 
in the bulge assuming a part icular mass-to-lig ht ratio 
(e.g. IMagorrian et al.l 119981 : [kormendv fc Richstond 
[l995t iGiiitekin et al.ll2n09b[r The M-a relation between 
SMBH mass and the velocity dispersion of the host 
galaxy also implies a physical coupling between forma- 
tion and growth of the black hole and i ts surrounding s 
(e.g iFerrarese & Mcrritt 2000; GcbhardEeFaD 120001 ). 
The driving mechanism behind these couplings and the 
M-a relation especially, is thought to be the result of 
mergers that drive accretion onto the SMBH, which 
can quench star formation as energ y released from 
the c entral engine drives the gas out (|Di Matteo et al.l 
[2005h . 

In particular, the study of an accretion disk around a 
SMBH begins with observations of the extended black- 
body spectrum emitted by the disk; a consequence of the 
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radial dependence of the temperature associated with the 
accretion disk. In SMBH accretion disks, this spectrum 
is thought to peak in the UV and is associated with the 
"Big Blue Bump" (jElvis et al.lll99l . Unfortunately, UV 
flux is extremely susceptible to scattering by dust and 
modeling to correct for this can induce large uncertain- 
ties. In X-rays, emission from inverse- Compton scatter- 
ing, magnetic flares and magnetic reconnection events 
associated with the accretion dis k are characterized we ll 
by a non-thermal power- law (e.g. iMcHardv et al.ir2004l l. 
Accordingly, X-ray flux can be another proxy for observ- 
ing accretion disks. 

In accreting systems, as material migrates toward the 
center, a fraction is also ejected into outflows that have 
both radiative and mechanical influences on their en- 
vironments. The exact physical nature has not yet 
been observationally determined, but outflows are seen 
in all types of accreting systems, i.e., proto-stellar ob- 
jects (e.g. iMundtl I1985D. neutron stars and stellar- 
mass black holes Ce.g. lMargon|[l98l . and SMBHs (e.g. 
iCohen et al.l Il979| ) . These outflows can reach super- 
sonic speeds when collimated into jets, eventually de- 
positing significant energy into their surroundings (e.g. 
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ICohen et al.l[l979l: iFabian et al.ll2002l: lAllen et al.1 12006') 

Material moving outward into their host galaxies also be- 
gins to cool via synchrotron radiation (e.g. Jones et al.l 
Il974f ). Observed in the radio frequencies, this non- 
thermal process emitted in the core of the system is pre- 
dicted to have a flat spectr um, independent of frequency 
(|Blandford fc Koiiigll Il979[ ) making it a great observa- 
tional tool for characterizing jet emission. 

Utiliz ing these two wave length regimes. iMerloni et al.l 
(■200a) . lFalcke et al.l (pOOl . and iGiiltekin et al.l (|2009al l 
have all suggested a "fundamental plane" of black hole 
activity connecting black hole mass. X-ray luminosity 
and radio luminosity. The plane spans over 9 orders of 
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magnitude in mass, 12 orders of magnitude in radio lu- 
minosity and 13 orde rs of magnitude in X-ray luminosity 
(|Merloni et al.|[2003l ) . This plane suggests that accretion 
(traced by X-ray luminosity) and jet production (traced 
by radio luminosity) are fundamentally linked together. 
Although the exact coupling is not understood, this re- 
lation implies accretion must be driving jet production. 

A known problem with the relation is the large scat- 
ter of the data abo ut the plane, (uradio = 0.88dex 
iGiiltekin et al.l[2009al ). This scatter can be attributed to 
observational errors or the result of non-simultaneouty 
between the observations themselves. Measuring the X- 
ray and radio luminosities at different times may have 
sampled different fluctuations in the accretion rate in 
individual sources, driving them away from the rela- 
ti on. The time between X-ray an d radio o bservat ions 
of iMerloni et al.l ([2003') and Giilte kin et al.l (|2009al ) are 
known to span a few years to a decade. 

To address these issues, and in order to examine disk- 
jet coupling at high mass accretion rates, we undertook 
a simultaneous X-ray and radio monitoring campaign of 
the Seyfert-1 AGN NGC 4051. This galaxy is relatively 
nearby (z=0. 002336), and the central black hole mass 
has been determined through re verberation mapping 
techniques ((1.91 ± 0.78) x lO^M© iPeterson et all [2001 
IDennev et ai][2009l (1.73 ± 0.55) x IO^Mq). NGC 4051 
is typically observed to accrete at approximately 5% of 
its Eddington luminosity (Peterson et al. 2004). The in- 
nermost orbital timescale of NGC 4051 is on the order of 
a few minutes to hours, defined as tdyn R/v^, where R 
is the radius assumed to be only a few gravitational radii 
from the black hole and is the orbital velocity. While 
the viscous timescales are on order of days to weeks for 
typical parameters, defined as t^s ~ tdynOt~^{H j R)^^ , 
where a is the vi scosity parameter in the s tandard a- 
disk prescription ([Shakura fc Sunvaevlll973() , and H is 
the scale-height of the disk. Variations in the accretion 
rate of NGC 4051 are thought to occur on these viscous 
timescales of a day to weeks. This is further supported 
by its highly variable spectrum, which is most notable in 
X-rays that vary up to a factor of 10 on weekly timescales 
(jUttiev et al.lll999D . The variability seen in NGC 4051 
was essential for our simultaneous monitoring campaign 
to probe different accretion rates. 

In this paper, we present the simultaneous X-ray and 
radio observations of NGC 4051, in effort to shed light 
on the implications of the fundamental plane of accretion 
onto black holes, and to explore jet-production in this 
Seyfert galaxy. 

2. DATA REDUCTION AND ANLAYSIS 
2.1. X-ray 

The Advanced CCD Image Spectrometer (ACIS) on 
the Chandra X-ray telescope was used to collect eight 
observations of NGC 4051 between 1 Jan 2009 (MJD 
54848.2) and 31 Jul 2009 (MJD 55043.1). We list these 
observations of approximately 10 ks exposures in Table [TJ 
X-ray observations are subject to photon pile-up. Pho- 
ton pile-up is when multiple low-energy photons arrive 
at the same CCD pixel in a single frame and register 
identically to one higher energy photon. It reduces the 
number of soft, i.e. low energy, photons that are detected 
while increasing the hard, i.e. high energy, photon counts 



which can mislead analysis by making the spectrum ap- 
pear harder than it actually is. To avoid photon pile-up, 
we used the mode with the minimal integration time, i.e. 
the Continuous Clocking mode, with integration times 
of only 2.85 ms. Such a fast integration time helps to 
prevent against multiple photons being counted as one 
higher energy photon, but does so at the expense of one 
spatial dimension. The resulting image is 1 x 1024 pixels 
with a resolution < 1", and an effective field of view of 
8.3' X 8.3', with spatial resolution in only one dimension. 

We used HEASDFT version 6.7, FTOOLS version 6.7, 
XSPEC version 12.5.1, CALDB version 4.1 and CIAO ver- 
sion 4.1 in the data reduction and analysis of these im- 
ages. The CIAD psextract routine was used to extract 
a spectrum and background for a point source as well as 
the companion auxiliary response functions (ARE) and 
response matrix functions (RME). To do so, we used 
a circular extraction region of diameter 2" centered at 
12h03m09.6s, -|-44°31'52'.'5, while another 2" diameter 
circular extraction region was used approximately 20" 
off-source for the background. This aperture size en- 
sured encapsulation of the source, assuming it was a 
point source and not elongated, but not so large as to 
include effects from the background. 

We analyzed X-ray data between the energy range of 
0.5-10.0 keV with a minimum of 10 counts per bin. A 
power-law component modeling the continuum was ini- 
tially fit to the data. To check for pile-up in our source 
before further modeling was attempted, we also used an 
annulus with inner radius of 0'.'9 and outer radius of 1'.'7 
to obtain an additional spectrum from each observation 
using psextract. Again, requiring a minimum of 10 
counts per bin, we fit power-law models to the annu- 
lar spectra for comparison. If pile-up was present, the 
spectra from the annuli would be softer than those ex- 
tracted from the central region. In general, we found 
that the two spectra were equivalent, within errors, sug- 
gesting that any pile-up is minimal in all observations. 

A description of the line-of-sight absorption is as fol- 
lows and was included in all modeling. The Galactic ab- 
sorption was modele d as an effective H c olumn density 
of 1.15x10^° cm~2 (jKalberla et al.l l2005| ) usii ig phabs 
Two absorption edges typical to Seyfert AGN () Reynolds! 
IT99I at K-shell rest energies of O VII and O VIII of 
0.739 and 0.871 keV respectively, were modeled using 
zedge. These three rest-frame values were frozen, while 
individual normalizations were allowed to vary. 

Eigure [1] shows a sample spectrum fit with a simple 
power-law model. This attempt at fitting the data re- 
sulted in a poor fit for all spectra and in particular the 
spectrum presented in Eigure[T]had a /v of 1770/366. 
By looking at the ratio of the data to the model, one 
can see the obv ious excess at low energies typical of 
Seyfert galaxies (|Revnold£lll997D . To better characterize 
the source flux, a disk blackbody component was added 
shown in Eigure [H Typically, it is thought that this soft 
excess is not the result of thermal blackbody component 
but of an a tomic process that is not tr ivial to model (e.g. 
iGierlinski fc Done 2004: iCrummv et a l. 2006). Eigure [3] 
plots the temperature and flux from the putative disk 
component. We find that the temperature does not vary 
with flux, suggesting it is not a blackbody thermal com- 
ponent. However, the addition of the disk blackbody 
does produce a formally better fit. 
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In Figure [21 we not only include a power-law compo- 
nent and a disk blackbody but also an unresolved Gaus- 
sian that models a narrow Fe Ka line and a broad Fe Ka 
line. The fit is improved to x^/j'=368. 3/353. It should 
be noted that the Fe Ka lines are not detected at more 
than the 2a level of confidence, except in the last two 
exposures. However, they were included in the modeling 
for completeness (see Table [2]). In Seyferts and AGN, 
Fe K-shell emission due to fluorescence and recombina- 
tion is the most prominent of the X-ray emission lines 
(e.g. [Miller 2007). The shape and broadening of the Fe 
Ka line, initially modeled by [Fabian et a l. ( 1989) for a 
zero-spin Schwarzschild black hole and [Laor, (.1991i) for 
a maximally spinning black hole, is dependent on the 
spin of the black hole, the inner and outer radius of the 
emitting region, the inclination angle of the disk and the 
efficiency of the disk emissivity. With this in mind, we 
modeled the broad Fe Ka line with XSPEC model laor 
that produces an emission line from an accretio n disk and 
includes general relativistic effects (|Laoiill99l[ ). Keeping 
the emissivity as a function of radius (oc R~^), the outer 
radius fixed at 400 GAI/c^, and the inclination fixed at 
30°, we let the inner radius vary. 

Finally, in order to compare to the fundamental 
plane of accre t ion on to black holes as described by 
iGiiltekin et"aL[ (|2009aD . just the power-law component 
within the energy range of 2-10 keV was used to com- 
pute the total flux associated with each observation. This 
component contributes to more than 95% of the total 
flux in the 2-lOkeV range, with the broad Fe Ka line 
contributing 2-5% of the total flux. We plot this light 
curve in Figure [4l and list the results for each spectrum 
in Table [T] The continuum flux varies by almost a fac- 
tor of 3, which can be attributed to fluctuations in mass 
accretion rate. 

2.2. Radio 

A combination of the Very Large Array (VLA) and Ex- 
tended Very Large Array (EVLA) antennae were used to 
observe NGC 4051 six times between 31 Dec 2008 (MJD 
54831.3) and 31 July 2009 (MJD 55043.1). These obser- 
vations were approximately 1 hour integrations centered 
at 8.4 GHz with a bandwidth of 50 MHz in two channels 
with two polarizations. We chose this frequency to en- 
sure self- absorption of the synchrotron emission was not 
a problem. Exposure dates and times are listed in Ta- 
ble [3] During each observation, the antennae switched 
from on source, located at lh203mls -|-44°31'8", to a 
phase calibrator located at 12h21m4s -1-44° 11'4" every 
3.33 s. 3C 286 was used as the flux calibrator with 3.33 
s of integration at the end of each observation. Dur- 
ing the seven months spanning the observations, the 
VLA/EVLA evolved from its A configuration, which has 
its longest baseline of 36.4 km, to its C configuration 
with a maximum baseline of only 3.4 km. This signif- 
icantly changed the resolution of the resulting images, 
and steps were taken to ensure a consistent comparison 
among the six observations. The resolution during each 
configuration is listed in Table [31 

The Common Astronomy Software Applications 
(CASA) package, version 3.0.0, developed by National Ra- 
dio Astronomy Observatory (NRAO) was used to reduce 
the radio observations. The routine setjy was used 



to set 3C 286 as the fiux calibrator, while the routine 
gaincal set the gain calibration incorporating both the 
phase calibrator source located at 12h21m4s -|-44°1'4" 
and 3C 286 as references. Either antenna VA06, VA08, 
VAIO, or VA12 was used as the reference antennae for 
these calibrations. Antennae that were off or showed 
discontinuities that were not characteristic of the overall 
data in either phase or amplitude were excluded, result- 
ing in only one or two EVLA antennae per data set. The 
CLEAN algorithm was then run to create an image, us- 
ing a threshold of 0.1 mJy as well as Briggs weighting 
and robust parameter of 0.5. Briggs weighting provides 
a smooth transition between natural and uniform visibil- 
ity weighting and is characterized by a robust parameter 
where 2.0 is approximately natural weighting and -2.0 is 
approximately uniform (Briggs 1995). No self-calibration 
was done because the source was too fain t, which is con- 
sistent with methods described by ,Giroletti fc Panessal 
{20M), who also observed NGC 4051 at 8.4 GHz. Fig- 
ure [5] shows an example of the type of images that were 
produced. 

After the CLEAN algorithm was run, the flux density 
was measured using CASA's imf it, which fit Gaussian 
curves to the peak intensities. Because of the four differ- 
ent configurations, we were only able to resolve structures 
in the A configuration, as seen in Figure[Sl which afforded 
the best resolution. In order to make a fair comparison 
between all the measurements, the Gaussian fits were re- 
stricted to an area of approximately 8" x 6", ensuring the 
inclusion of all the structures presen t in the images. For- 
tunate ly, as seen in the work done bv lGiroletti fc Panessal 
(|2009() . there does not appear to be any extensive large 
scale diffuse emission that would be resolved out in higher 
resolution configurations. If this were the case, the flux 
would decrease with higher resolution, when in fact Ta- 
ble [H demonstrates that the flux in the B configuration, a 
higher resolution, is larger than the fiux in the C config- 
uration, a lower resolution. This strongly suggests that 
the radio variability observed between the B and C con- 
figurations is intrinsic, not instrumental. However, it is 
unclear if the diffuse emission has been resolved out in 
the A configuration. For this reason, we quote our results 
with and without the A array measurements. The mea- 
surements including the integral fiux densities as well as 
the peak fiux densities are listed in Table [3] and the light 
curve is shown in Figure [6l The errors in these measure- 
ments also include a systematic error of 3% of the total 
fiux, added in quadrature, to account for the calibration 
errors. These measurements do vary by a factor of 3, 
which is co mparable to the X-ray v ariations. In order to 
compare to IGiiltekin et al.l (|2009aD , we scaled the fluxes 
as oc from 8.4 GHz to 5 GHz as|Ha (2002) did. 

Inspection of the image produced in the A configu- 
ration reveals evidence of jet production. As shown in 
Figure [3 three structures are present: a central source 
and two extended lobes. The central radio lobe is likely 
associated with the black hole, while the northwest and 
southeast lobes may be associated with the endpoints of 
the outfiows in this system. The alignment and vary- 
ing brightness between the three lobes suggests this is 
a jet structure where the jet is partially projected into 
our line-of-sight, placing the northwest lobe closer and 
less obscured then the southeast lobe. We also note 
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that these lobes are not seen in the antennae psf and are 
therefor e not artifacts of the Clean algorithm. Observa- 
tions b vlChristopo ulou et aI](ll997D . iGiroletti fc Panes"sal 
(|2009h . and Kukul aeTalJ (|1995D note similar structures 
as well as the potential evidence for jet production. 

2.3. Comparison 

These measurements by themselves can give insights 
into the physical processes that occur in AGN, but 
through a comparison of the X-ray and radio luminosities 
we can explore how accretion is coupled to jet produc- 
tion. To quantify any correlation between the X-ray and 
radio fluxes, we calculated the Spearman's rank coeffi- 
cient (p). This p describes how well a data set exhibits 
a monotonic behavior between two variables; p < cor- 
responds to variables that are anti-correlated. We find a 
p of —0.66, with a probability of only 16% of zero corre- 
lation. 

Figure [7] plots the X-ray vs. radio luminosities. We as- 
sumed a distance of 10.0 Mp c derived from the redshift, 
z = 0.002336 (P eterson et a l. 2004), and cosmological 
parameter Ho— 70 km Mpc to obtain these luminosi- 
ties. Figure [7] also includes the least-squares fit of a first 
order polynomial to the data points, described by the 
following relation: 

log (-0.72±0.04)logLx-ray + (64±2) (1) 

To test the robustness of the anti-correlation, we used an 
F-test to compare the results of the first order polynomial 
given above to fits with a zeroth order polynomial as well 
as the fundamental plane relation of iradio oc Lx-tay • 
At an 81% confidence level, we were able to rule out 
the flat model. However, at a 96% confidence level, the 
fundamental plane relation is excluded by our model. 

We also looked at the correlation between the data 
points excluding the A configuration, in the case that 
this array resolved out a substantial amount of flux. 
We found that the correlation was still inversely propor- 
tional, described as, 

logiradio = (-0.12 ± 0.05) logix-ray + (40 ± 2) (2) 

Using the same F-test, our second model may be consis- 
tent with being flat, for it only excludes a flat distribution 
at the 60% confidence level. However, the positive fun- 
damental plane relation was again excluded, this time at 
a 98% confidence level. 

To place the data from NGC 4051 in context of the fun- 
damental plane of accretion onto black holes, we plotted 
our data and be st fit line against the da ta and the best fit 
line described in lGiiltekin et al.l (l20q9aD in F igurelSl The 
best fit line found by lGiiltekin et al.l ()2009al ) is described 
as follows. 



logiradio = (0.67 ± 0.12) l0gLx-ray + 

(0.78 ± 0.27) logMBH + (4.80 ± 0.24) 



(3) 



One should notice that our data are not described 
by the positive correlation of Lradio oc ixSy 
noted by the solid line, but instead by the dashed line 
and dotted line described by Lradio fx ix-ray'^ ^^^'^ 



o.i2±o.05^ respectively. In fact, our relation 
differs from this relation by 11.6ct, while including the 



-t'radio C< ^X-ray 



A configuration, and 6.5(7, while not including it, when 
dividing the difference between the power-laws by the 
larger of the two errors. The fundamental plane is known 
to have a large scatter, and NGC 4051 lies close to the 
plane. This exercise, however, makes clear how different 
the X-ray and radio coupling in NGG 4051 appears to 
be. 

Finally, we attempt to quantify whether or not there is 
a significant time delay between fluctuations in the X-ray 
and radio luminosities, which can be seen in Figures |4]& 
[6] respectively. We used a discrete cross-correlation func- 
tion (DCCF) as described by Edel son fc KrolikI ()1988[ ) 
to quantify this delay. The sparse and uneven sam- 
pling of the all six observations poses a natural limi- 
tation on the results of this analysis. This led us to 
linearly interpolate the data within the seven month pe- 
riod with uniform spacing. The time steps between in- 
terpolated points were 17 days, cor responding to the 
shorte st tin ie between observat i ons. I White fc PetersonI 
([1991 and IGaskell fc PetersonI (flgSTT describe similar 
techniques for interpolating and cross-correlating data in 
their respective variability studies of AGN. The results, 
shown in Figure HI place the minimum of the DCCF of 
-0.48 ±0.1 at -2.5±5.3 days, suggesting that X-ray dips 
are leading the radio flares. 

3. DISCUSSION 

In this paper, we present eight Chandra X-ray observa- 
tions and six VLA/EVLA radio observations of Seyfert- 
1 NGC 4051 taken over a seven month period. By si- 
multaneously measuring in X-ray and radio bands ev- 
ery two to four weeks apart, we are able to probe vari- 
ations in the accretion rate that occur on the viscous 
timescale of a few days to weeks in NGC 4051. The 
observations reveal significant variability in both X-ray 
and radio bands. The first 8.4 GHz observation on MJD 
54831.3 also shows evidence of jet production in the form 
of two distinct radio lobes, contrary to the idea that ra- 
dio quiet galaxies like NGC 4 051 lack jet production. 
In fac t, work by iFalcke et al.l ([2001,) and Nagar et aX] 
()2002f ) suggest that elongated, non-thermal emission is 
common in Seyfert and LLAGN. The lobes appear to be 
resolve d into compact impact reg ions at very high reso- 
lution (jGiroletti fc Panessall2009l ). A variability analysis 
shows that the 2-10 keV X-ray luminosity and 8.4 GHz 
radio luminosity are inversely correlated according to 
iradio oc ivx° Jfy=^° °^ This diffcrs by 11.6cr from the cur- 
rent fundamental plane relations which correlates the ra- 
dio luminosity to X-ray luminosity as Lradio « ^xSy 
for a fixed mass ()Giiltekin et al.l l2009aD . Furthermore, 
if extended emission is resolved out in the highest res- 
olution image, then excluding the A configuration data 
point, the correlation is iradio « ^x-ray'^ This still 
differs from the fundamental plane by more than 6a and 
is not consistent with the relation. 

The results help to shed light on the fundamental 
plane of accretion onto black holes fMerl oni et al.l [20031 : 
IFalcke et al] 120041 : LGiiltckin et al. 2009a,). This plane 
demonstrates that accretion is linked to jet production, 
as evidenced by the positive correlation between X-ray 
and radio luminosity. The measurements in the funda- 
mental plane are generally not simultaneous, and one 
goal of this study was to understand if simultaneity be- 



5 



tween X-ray and radio measurements would reduce the 
scatter in the plane. NGC 4051 was found to lie near 
to the fundamental plane, but showed a negative corre- 
lation between X-ray and radio luminosities, contrary to 
the positive correlation of the plane. The inverse correla- 
tion suggests that some systems may lie near to the fun- 
damental plane, but vary like NGC 4051, moving across 
the plane and thus increasing the total scatter. 

At least a few other case studies have also shown a 
separate and inverse correlation deviating from the posi- 
tive correlation suggested by the fundamental plane. The 
first is 3C 120, a SMBH with a mass of approximately 
5.5 X IO^Mq that is observed to acc rete at approximately 
10% LEdd ()Peterson et all I2004D . iChatteriee et al.1 
(|2009D present a five year study of 3C 120 using RXTE 
and VLBI to obtain X-ray and radio data respectively. 
They use a discrete cross-correlation function to describe 
the correlation between the X-ray at 2.4-10 keV and the 
radio at 37 GHz. They find that the greatest amplitude 
in the DCCF to be —0.68 ±0.11, corresponding to an in- 
verse correlation at a 90% confidence level. This ampli- 
tude is consistent with the results of the cross-correlation 
analysis in NG C 4051. The minimum i n of the DCCF of 
3C 120 cited in lChatterjee et all (|2009| ) corresponds to a 
time lag of 120 days ±30 days, with the X-ray dips lead- 
ing leading the radio flares. This time delay corresponds 
to approximately 4 days in NGC 4051 when scaled us- 
ing their respective masses, which is consistent with our 
data. 

At the opposite en d of the mass scale, the 14M0 black 
hole GRS 1915±105 ()Greiner et al.ll2O0l[ ) sometimes also 
shows an inverse relatio n between simu l taneo us X-ray 
and radio observations. iRau fc Greineii (|2003[ ) present 
a survey of X-ray and radio observations using RXTE 
and the Ryle Telescope between 1996 November to 2000 
September. The 1-200 keV X-ray fiux showed no corre- 
lation to the 15 GHz radio observations. However, the 
20-200 keV continuum did show an inverse correlation to 
the 15 GHz fiux described by a Spearman's rank-order 
correlation coefficient of -0.75. The hard X-ray band in 
GRS 1915+105 excludes direct emission from the disk, 
which is consistent with using X-ray emission instead UV 
emission in NGC 4051 and 3C 120. 

The prototypical stellar-mass black hole, Cygnus X- 
1, may also show an inverse trend at high luminosity. 
iGallo et al.l (|2003[ ) study the disk-jet connections in this 
stellar-mass black hole by analyzing RXTE All Sky Mon- 
itoring X-ray data from 2-11 keV and the Ryle Telescope 
radio data at 15 GHz between 1996 January and January 
2003. Cygnus X-1 does follow LRadio oc L^Iray until 
it reaches approximately 2% of its Eddin gton lum i nosity 
when the radio flux density turnovers. G allo et al.l (l2003f ) 
describe departures from this relation at high X-ray lu- 
minosity as quenching of its jet production as Cygnus 
X-1 moves into its "high/soft" state. 

Observations of NGC 4051, 3C 120, GRS 1915+105, 
and Cygnus X-1 all show that X-ray and radio flux fol- 
low an inverse relation when observing at nearly simul- 
taneous times when the sources are emitting at 1-10% 
of Eddington. Given that an inverse correlation is seen 
in a quasar, a Seyfert, and two stellar-mass black holes 
at high Eddington fractions, it is possible that a dis- 
tinct mode of disk-jet coupling holds at high Eddington 



fractions. This would go beyon d a simple quenc h ing o f 
jet production, as discussed bv iMaccarone et all (|2003[ ) 
and lGallo et all ([200l , since the jet production does not 
turn off entirely (evidenced by continuous radio emission 
and jet structur es, especially in NG C 4051 and 3 C 120; 
iGiroletti fc Panessa .2009,: Chatteriee et al.l[2009( ). This 
is the first study to probe a Seyfert galaxy in X-ray and 
radio bands on viscous timescales of its inner disk. In the 
future, we will undertake an extended monitoring cam- 
paign of NGC 4051 to further characterize this relation 
as well as determine any time lags between X-rays and 
radio fluxes. 

The authors gratefully acknowledge Joan Wrobel, Phil 
Uttley, and Ian McHardy, and the anonymous referee 
for their insights and comments that improved this pa- 
per. They gratefully acknowledge Steven T. Myers for 
his help with the CASA software package. J. M. M. grate- 
fully acknowledges support through the Chandra guest 
observer program. S. M. gratefully acknowledges support 
from a Netherlands Organization for Scientific Research 
(NWO) Vidi Fellowship and from the European Commu- 
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under grant agreement number ITN 215212 "Black Hole 
Universe". E. M. C. gratefully acknowledges support 
provided by NASA through the Chandra Fellowship Pro- 
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NOTE ADDED IN PROOF: While our paper was be- 
ing reviewed, a separate paper on radio and X-ray ob- 
servations of NGC 4051 was accepted for publication 
i jJones et al.l I2010D . Our paper was accepte d for pub- 
lication only days later. iJones et al.l (|20100 report on 
many more radio observations, and employed RXTE to 
obtain X-ray flux points. Treating th e possibility of dif - 
fuse nuclear emission with great care, IJones et al.l (|20Tof ) 
find a nearly fiat radio-X-ray correlation. Our results are 
broadly consistent with that finding when the EVLA ob- 
servation in the A configuration is excluded. The main 
advantages of our work are the close timing of radio and 
X-ray observations, and the ability to separate distinct 
X-ray flux components in Chandra spectra. Both papers 
show that NGC 4051 lies significantly below the Funda- 
mental Plane. Within the context of other black holes 
accreting at high Eddington fractions (a point of em- 
phasis in this paper), both sets of results support the 
possibility that the coupling between the disk and jet in 
NGC 4051 may be different than in low- luminosity AGN 
and other extremely sub-Eddington sources. We thank 
Ian McHardy and Phil Uttley for helpful conversations 
regarding NGC 4051 and difi^use nuclear radio flux. 
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TABLE 1 
X-RAY Observations 



Date of 


Exposure 


Count 


r 


r Flux 


Disk Blackbody 


Disk Blackbody 




Observation 


Time 


Rate 




(z-lUkeV ) 


Temperature 


llux (2— lUkeVj 




(MJD) 


(ks) 


(counts s^-"-) 




(10^"'""'" erg s~i cm^'^) 


(keV) 


(10^-'^^ ergs s~l cm~-^) 




54838.2 


10.2 


15.6 ± 0.04 


r, r,-+0.03 

^•^^-0.03 


n 97+0.09 
-0.09 


n 10+0.02 
'^•^''-0.02 


10.7+^-? 


646.1/426 


* 54874.9 


1.1 


16.3 ± 0.12 


2 39+0-11 


9 co+0.17 
Z.0c5_g 20 


17+0-04 

^■^'-0.03 




244.6/213 


54898.3 


10.1 


5.3 ± 0.02 


-I 7O+0-08 
^•'■'-0.03 




17+0-01 
'^•^'-0.02 




394.0/381 


* 54943.1 


10.1 


4.7 ± 0.02 


2 

^■I'^-O.OT 


1 14+0-05 
^■^^-0.08 


n 10+0. 03 

U-lS-0.03 


7+4.4 


368.0/353 


54988.0 


10.1 


4.6 ± 0.02 


1 OO+0-06 
^■'-''-'-0.07 


^■°'-0.12 


19+0-01 
'-'•1^-0.01 


16.0 til 


370.0/394 


55005.8 


10.1 


5.6 ± 0.02 




-I 90+0.07 

^■^'-'-0.08 


n is+o-oi 


S 1+3-2 


501.8/364 


55025.4 


10.1 


13.1 ± 0.04 


r, 9Q + 0.03 


9 C.7+O.O6 
-0.07 


n 10+0.02 


5 1+4-1 


578.7/411 


55043.1 


10.1 


14.9 ± 0.04 


9 9Q+0.03 


9 07+0-07 
-0.06 


17+0-02 
'-'•l'-0.02 




627.7/420 



Note. — This table shows the X-ray observations made by Chandra in the continuous clocking mode. The data was modeled using 
Xspec and unabsorbed fluxes for both the power-law and disk blackbody are presented here. The power-law used in the analysis varies by 
a factor of 3. The * refers to the observations that were not used in this analysis because of the lack of simultaneous radio measurements. 



TABLE 2 

X-RAY Fe Ka LINE 



Date of 
Observation (MJD) 



Narrow 
Fe Ka Equivalent 
(keV) Width (keV) (10" 



Flux 



Broad Fe Ka 
line (keV) 



Rin 
(GM c- 



Broad 
Equivalent 
Width (keV) 



(10 



-13 



Flux 
erg s~^ cm" 



54838.2 

* 54874.9 
54898.3 

* 54943.1 
54988.0 
55005.8 
55025.4 
55043.1 



1.40 



-0.04 
-JVA 



0.04: 



-0.04 
-0.04 



1.36 



+1.29 
1.26 



4P.+0.09 
■^^-0.07 


04+°-°'' 





84 


zin+0-05 

■40-0.04 


14+0-07 
'^•1^-0.07 


1 


81 


Qn + 0.05 

.39_^^ 


06+° °^ 

"•'-'"-0.05 


1 


32 


QR + 0.03 




1 


73 


zin+0-11 




1 


02 


(,7+0.07 
■t''-0.08 


0.02tH? 





64 



+ 1.31 
-NA 
+0.96 
-0.98 
+ 1.20 
-1.25 
+ 1.05 
-1.44 
+ 1.31 
-NA 
+ 1.45 
-0.15 



6 95+^^ 
O-yO_0.61 

6.34+0/i 
6.62l^3l, 

6-9llS'5li 
6 35+0-13 

6 43+0-13 
6 34+0-15 



K+20 
-JVA 
-165 
-JVA 
hl3 
-JVA 
1 9+69 
^■^-NA 
3+31 

-JVA 
1+20 
-32 



5.6^ 
25^ 

3.r 



3.31 
44! 
20^ 



31+0-24 
"-•^1-0.18 
54+0-80 

15+0-38 

7n+o « 

"■ '^^-0.45 
39+0.26 
U.03_Q 28 
o O9+0.26 
"-'-'^-0.27 

0-30t°;l? 

0.29«;}^ 



7+5.48 
-4.11 
-15.0 
-9.18 
7+5.76 
-JVA 
7+3.53 
-3.66 

6 59+*-^^ 

U->J»_4.75 
1+3.20 
-3.23 
-1+2.72 
-2.49 
-,+2.98 
-2.86 



7.07Z 
10.1^ 

2.27^ 
5.67^ 



3.94: 
6.6OI 
7.29^ 



Note. — This table shows the two Fe Ka lines included in the X-ray models. To calculate the narrow Fe Ka lines, zgauss with energy 
restricted between 6.35-6.97, a width of a =0 (unresolved), redshift of z=0. 002336 and varying normalization was used. For the broad Fe 
Ka line, laor the energy was restricted to vary from 6.34-6.95 keV, the power-law dependence was frozen at 3, the inner edge was allowed 
to vary from 1.235 to 400 GM/c^, the outer radius was frozen at 400 GM/c^, the inclination angle was frozen at 30° and the normalization 
was allowed to vary. All the lines were at least marginally detected and therefore included in modeling techniques, except the narrow Fe 
Ka on MJD 54874.9, where it registered flux for the best fit parameters. The * refers to the observations that were not used in this 
analysis because of the lack of simultaneous radio measurements. 



TABLE 3 
Radio Observations 



Date of Observation 


Exposure Time 


Configuration 


Resolution 


Flux Density 


Peak Flux Density 


(MJD) 


(ksec) 




(arcsec) 


(pc) 


(mjy) 


(mJy) 


54831.3 


3.6 


A 


0.3 X 0.2 


15 X 10 


1.73 ± 0.06 


0.65 ±0.02 


54899.2 


3.5 


B 


0.9 X 0.7 


44 X 34 


5.99 ± 0.20 


1.13 ±0.04 


54987.1 


3.6 


BnC 


2.2 X 0.9 


110 X 44 


5.66 ± 0.20 


1.47 ±0.05 


55005.2 


3.6 


C 


3.1 X 2.3 


150 X 110 


4.97 ± 0.17 


2.22 ±0.08 


55027.0 


2.4 


C 


2.5 X 2.0 


120 X 97 


4.99 ± 0.20 


2.08 ±0.08 


55043.1 


3.6 


c 


3.1 X 2.2 


150 X 110 


4.78 ± 0.16 


2.22 ±0.08 



Note. — This table gives the radio observations taken at 8.4 GHz with a 50 MHz bandwidth in two channels. The VLA/EVLA was in 
4 different configurations during the entire campaign, s tart i ng in the A with longest baseline and ending at C with the shortest baseline. 
In the analysis, we scaled the observations to 5 GHz as lHol 120021) did using the total flux densities given in column 5. 
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Energy (keV) 

Fig. 1. — This is a sample spectrum from MJD 54943.1 modeled with just a power-law component, in red, as the Xspec model 
phabs (po) *zedge*zedge. The Galactic absorption was modeled as an effective H column density of 1.15x10^'^ cm~^ l IKalberla et al.l 
1200511 and the K-shell absorption edges of OVII and OVIII were frozen at their respective rest energies of 0.739 keV and 0.871 keV. We 
initially fit the spectrum from 2-10 kcV and then extended the model to lower energies. The fit produced a x^/i'= 1770/366. 




Fig. 2. — This is the same spectrum as Figure [T] but now includes a power-law, disk blackbody, a narrow Fe Ko line and a broad Fe 
Ka line, modeled with Xspec model phabs (po+diskbb+zgauss+laor)*zedge*zedge. We used the same Galactic absorption and absorption 
edges as in Figure [T] The power-law is the red dashed line. The disk blackbody is the blue dot-dashed line. The narrow Fe Ka line is the 
green dotted line. The broad Fe Ka line is the black triple dotted-dashed line. Finally, the solid black line is the sum of all the components. 
The fit produced a x^/i'=368.0/353. 



8 





' 1 


i 1 1 1 1 1 1 1 i 1 1 1 1 1 
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1 


1 .... 1 ... . 



Diskbb Flux (ergs s"' cm"^ ) 

Fig. 3. — This figure siiows tlie comparison between the soft-excess disk blackbody temperature and its blackbody flux. A constant 
temperature independent of flux as seen here suggests a purely phenomenological interpretation of the disk blackbody component. 
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5.48'10^ 5.49'10^ 5.50«10^ 5.51 '10^ 
Modified Julian Date 



Fig. 4. — This is the X-ray light curve taken by Chandra in the continuous clocking mode between 2-10 keV. It shows a variability of a 
factor of 3. The points in red do not have simultaneous radio observations and are not used in the correlation analysis. 
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J2000 Right Ascension 

Fig. 5. — Shown here is a radio observation made by VLA/EVLA in the A-configuration on 31 Dec 2009 (MJD 54831.3). The contours 
are [0.2, 0.4, 0.6, 0.8] x 0.6 mjy beam~^. The beam pattern (3" x 2") is also shown in the lower left corner. The center is associated with 
the black hole, while the extended lobes northwest and southeast are associated with the endpoints of jets of the system. 
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Fig. 6.— 



This figure shows the 8.4 GHz light curve taken by the VLA/EVLA. The radio flux density varies by more than a factor of 3. 
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35.5 




41.0 



41.2 



41.4 41.6 41.8 

log(Lv_R„„ [ergs s 'l) 

Fig. 7. — This figure shows the inverse correlation between ??-ray^af?3JTaaio mminosixies. We used the continuum flux from 2- 
and we scaled the radio flux as Fy oc u~^'^ (Ho 2002) from 8.4 GHz to 5 GHz in order to correctly compare to [Giiltckin ct al. 
The dashed line is a fit to our data given by the relation, logL^adio = (—0.72 ± 0.04) logLx-ray + (64 ± 2). The dotted line excludes the 
A configuration data point, and is described by the relation, logL^adio = (—0.12 ± 0.05) logLx-ray + (40 ± 2) . 



-10 keV, 
(2009a). 



11 



42.5 



39.7 - 



CO 

CD 



" 36.9 - 



o 
EE 



D) 
O 



34.0 - 



31.2 



LINERS/LLAGN 
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31.5 



33.8 
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0.67log(Lx_Rav [ergs s''l)+0.78log(M [Mcoi-. . 

Fig. 8. — This is a projection of the fundamentarplaire Br%.ccrotion ontoDiack holes as acscriped pv ^ filtckin et ahl H2009al ) ). The red 
points are Seyfert galaxies, while blue points are LINERS or LLAGN. The solid black line is Giiltekin et aU I l2009al) best fit fine. Our data 
points are shown in black, with our best fit as the dashed line and dotted lines as described in Figure [T] 
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Fig. 9. — This figure is of the discrete cross-correlation function (DCCF) between X-ray and radio luminosities. We fit the DCCF with 
a quadratic function, and the minimum occurs at -2.5±5.3 days with a DCCF value of -0.48±0.1. A negative delay implies that the radio 
is lagging the X-ray. 
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